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ABSTRACT 

We investigate accretion disc models for the X-ray emission of Seyfert-1 galaxies and 
the hard state of black-hole binaries. We concentrate on two hot accretion disc models: 
advection-dominated accretion flow (ADAF) and recently found luminous hot accre- 
tion flow (LHAF). We solve for the global solution of both ADAF and LHAF to obtain 
the electron temperature, T e , and Thompson optical depth, r, at the radius where most 
of the radiation comes from. We adopt two kinds of electron energy equations. In one, 
only synchrotron and bremsstrahlung radiation and their Comptonization are consid- 
ered. The other is parameterized by a constant Compton parameter to model the case 
in which thermal Comptonization of external soft photon is important. We compare 
the calculated T c and r with the observational values obtained by fitting the average 
spectra of Seyfert-1 galaxies and black-hole binaries using thermal Comptonization 
model. We find that the most favoured model is an LHAF with parameterized elec- 
tron energy equation, with ADAFs predicting too high T c . Also, the LHAF, but not 
ADAF, can explain large luminosities in excess of 10 per cent of the Eddington lumi- 
nosity seen in the hard state of transient black-hole binaries. 

Key words: accretion, accretion discs - black hole physics - galaxies: active - galax- 
ies: nuclei - hydrodynamics - radiation mechanisms: thermal 



1 INTRODUCTION 

X-ray spectra of black-hole X-ray binaries in the hard state 
and Seyfert-1 galaxies usually consist of a power-law with 
a high-energy cutoff, Compton reflection, and Fe Kq emis- 
sion (e.g., Nandra & Pounds 1994; Zdziarski et al. 1995; 
Nandra et al. 1997). The cut-off power law is produced, 
most likely, by thermal Comptonization in a hot, mildly 
relativistic, plasma (see, e.g., Zdziarski & Gierlinski 2004 
for a review). Physical models for the hot plasma include a 
hot accretion disc, which we will discuss below, a magnetic- 
dominated corona (e.g., Galeev, Rosner & Vaiana 1979; 
Haardt & Maraschi 1993; Beloborodov 1999a), and an ion- 
illuminated accretion disc (Spruit & Haardt 2000; Deufel & 
Spruit 2000). Among them, the hot accretion disc model has 
the most clear dynamics and fewest free parameters, and we 
concentrate on that model in this paper. 

Shapiro, Lightman & Eardley (1976, hereafter SLE) 
proposed the first hot accretion-disc solution. The available 
gravitational energy is converted in the process of viscous 
dissipation into the thermal energy of ions, which are much 
heavier than electrons. Since the only coupling between ions 
and electrons is Coulomb collisions which is rather weak, 



and the radiation of electrons is much stronger than that 
of ions, the ion temperature, Ti, is much higher than that 
of electrons, T c . The Thomson optical depth of the flow, r, 
is low, so T e can be as high as 10 9 K, high enough to pro- 
duce X-ray photons. The SLE solution is thermally unstable 
(Pringle 1976). In addition, energy advection, which can be 
very important, is neglected in the energy equation of ions 
in SLE. 

The second hot accretion-disc solution, advection- 
dominated accretion flow (ADAF; Ichimaru 1977; Rees et 
al. 1982; Narayan & Yi 1994, 1995; Abramowicz et al. 1995; 
see reviews by Narayan, Mahadevan & Quataert 1998; Kato, 
Fukue & Mineshige 1998), does include advection, which 
dominates the energy transfer for ions. The ion energy equa- 
tion reads g adv = q v i B — qi e , with g adv , 9vis and qi e being 
the rates of energy advection, viscous heating and Coulomb 
cooling, respectively. In a typical ADAF, due to t < 1, 
qic -C 9vis ~ i/adv, fe., the viscous heating is balanced by 
advection cooling. We define a dimensionless accretion rate 
as 
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where M is the mass accretion rate, and Le is the Eddington 
luminosity (note that in this definition we have not included 
the factor of 0.1 used by, e.g., Narayan & Yi 1995). With 
the increase of rh, since gi c oc rh 2 whereas g v j a oc rh, Qj e in- 
creases faster than q v i s , Coulomb cooling becomes more and 
more important. When rh reaches a critical value, denoted 
as rhi here, we will have qvis ~ qic, so advection fails to 
be the dominant cooling mechanism. We call rhi the criti- 
cal rate of ADAF (Narayan, Mahadevan, & Quataert 1998). 
The existence of rhi and their low radiation efficiency make 
ADAFs rather dim. Therefore, this solution is very success- 
ful in explaining low luminosity /quiescent states of black- 
hole binaries and low-luminosity AGNs. However, it is not 
clear whether it applies to luminous X-ray sources such as 
Seyfert-1 galaxies and the luminous hard state of black- hole 
binaries. 

Then, Yuan (2001, hereafter Paper I) found a new hot 
accretion solution above mi, the so-called luminous hot ac- 
cretion flow (LHAF). We emphasize that the equations de- 
scribing LHAF are completely identical to those of ADAF, 
so LHAF is actually a natural extension of ADAF to ac- 
cretion rates above rhi. In an LHAF, rh > rhi, so we have 
q v is < qic- The reason why hot solutions still exist above rhi 
is that the compression work, q CO m, provides additional heat- 
ing in addition to q v i B , so q coul + q v i B > q\ c (see Paper I for 
details) if rh is not too large. Since g a d v = q v is — qic < 0, ad- 
vection plays a heating rather than a cooling role in LHAF. 
In other words, the entropy of the accretion flow in an LHAF 
is converted into radiation, similar to the cases of spherical 
accretion and cooling flow in galactic clusters 1 . 

A unified description of the three hot accretion disc 
models, namely the SLE solution, ADAF and LHAF, has 
been presented in Yuan (2003). Yuan (2003) has also shown 
that if thermal conduction is neglected, LHAF is thermally 
unstable under local perturbations. However, the timescale 
of the growth of the thermal perturbations is in general 
longer than the accretion timescale if rh is not too large, 
despite the relatively strong radiative cooling in LHAF, so 
the solution can survive. Furthermore, inclusion of thermal 
conduction generally has a very strong stabilizing effect. 

From ADAF to LHAF, both the accretion rate and the 
radiation efficiency increase continuously (but the efficiency 
is still lower than in the standard thin disc), so LHAFs are 
significantly more luminous than ADAFs. Thus, LHAFs of- 
fer a promising solution to the problem of the nature of 
X-ray emission of luminous accreting black holes. Here, by 
comparing the predictions of the ADAF and LHAF solu- 
tions for the T e , r, and the bolometric flow luminosity, L, 
with results of fits to hard X-ray observations, we investigate 



1 We expect the existence of solutions which are ADAFs far from 
the black hole and become LHAFs at smaller radii when rh is not 
very large, because the radiative cooling rate, q ra< j, is very low 
at large radii. Such solutions were not presented in Paper I. We 
find that they do exist, but only in a narrow range of rh. This 
narrowness is, in fact, due to the definition of LHAF as satis- 
fying g v ; s < <jfj c . If, instead, we define LHAF solutions as those 
satisfying q v - ls < q ra d> we would have found the corresponding 
ADAF-LHAF transition solutions in a wider range of rh. This 
is because at R> 10 2 R S , q is 2> q r ad> where R s = 2GM/c 2 . On 
the other hand, the two cases become identical at small radii, 
R<10 2 R S . 



which of these two solutions is more appropriate to describe 
the nature of X-ray emission of luminous black hole sources. 

X-ray emission of Seyfert galaxies and black-hole bina- 
ries in the hard state is most likely due to thermal Comp- 
tonization of soft photon by thermal hot electrons. For a 
given geometry, the thermal-Comptonization spectrum is 
determined by two parameters of the hot plasma, the elec- 
tron temperature, T e , and the Thompson optical depth, r. 
Fits with this model, summarized in Section |5] below, pro- 
vide observational constraints on T c and r, as well as on 
L. We then compare them with the values predicted by the 
ADAF and LHAF global solutions. 

On the theoretical side, the exact cooling mechanism 
of electrons remains unknown. In particular, Comptoniza- 
tion of synchrotron radiation photons may be not the dom- 
inant process in luminous X-ray sources. Evidence support- 
ing this point includes a correlation between the strength of 
Compton reflection and the X-ray spectral index (Zdziarski, 
Lubihski, & Smith 1999; Zdziarski et al. 2003) and theoreti- 
cal arguments showing that this process is often not capable 
to provide the required flux of seed photons required (SLE; 
Zdziarski et al. 1998; Wardzihski & Zdziarski 2000). 

Sources of soft photons other than the synchrotron 
emission are, however, very likely to exist. The soft photons 
can come from a cold disc just outside of, or penetrating 
into, an inner hot accretion flow, if the truncation radius 
of the cold disc is not too large (Poutanen, Krolik & Ryde 
1997; Zdziarski et al. 1999). Another possibility is that the 
accretion flow consists of two phases, with some cold clouds 
or clumps embedded in the hot accretion gas (Guilbert & 
Rees 1988; Celotti, Fabian & Rees 1992; Kuncic, Celloti & 
Rees 1997; Krolik 1998) . In this case, the emission from the 
clumps can serve as the seed photons of Comptonization. 
There is also very strong observational evidence for cold 
medium extending close to the black hole in Seyferts and 
black-hole binaries from relativistic broadening of the reflec- 
tion/reprocessing features (e.g., Nandra et al. 1997; Zycki, 
Done & Smith 1998). In accord with those observations, the 
LHAF model predicts that when rh is very high, the hot 
flow can collapse and form an optically-thick disc in the in- 
nermost region (Paper I), which disc will also be a copious 
source of soft photons. 

For those reasons, the determination of the electron en- 
ergy equation can be very complicated. An approach to over- 
come this difficulty is to parameterize the energy equation 
of electrons using the Compton parameter, 

y = 4r^, (2) 

where r = <jTn c H corresponds to the disc scale height, H, 
<tt is the Thomson cross section, m c c 2 is the electron rest 
energy, and n c is the electron density. A given value of y 
corresponds quite well to the X-ray photon index V (Ghis- 
ellini & Haardt 1994; Poutanen 1998; Beloborodov 1999b), 
which is observationally well constrained. Note that some 
other definitions of y also exist, modifying the above ex- 
pression for the cases of large r and relativistic T e . In par- 
ticular, Beloborodov (1999b) used the definition of y = 
4[kT c /m c c 2 + 4(kT e /m e c 2 ) 2 ](r + T 2 ), for which case he found 

I-^r 2/9 (3) 

in a spherical geometry. 
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Furthermore, y (or, equivalently, V) is often almost con- 
stant with varying flux for luminous black-hole sources (e.g., 
Cyg X-l, Gierliriski et al. 1997; GX 339-4, Zdziarski et al. 
1998, 2004; Wardzihski et al. 2002; IC 4329A, Fiore et al. 
1992), which tells us that the relative cooling rate is typi- 
cally more stable than the X-ray flux from the source. Then, 
instead of considering the various sources of soft photons, 
we replace the electron energy equation by assuming that 
there is a local source of seed photons leading to a given 
parameter y, and solve for the disc equations. This is in fact 
the approach adopted by SLE (who assumed y = 1), and by 
Zdziarski (1998), who generalized the hot disc model of SLE 
by including energy advection in the ions energy equation 
and allowing y to be a free parameter. 

We present the data we use in Section |5] In Section |3] 
we solve the radiation-hydrodynamic accretion equation pa- 
rameterized by Compton parameter, y. For completeness, we 
also investigate the standard case, where only Comptoniza- 
tion of synchrotron and bremsstrahlung photons is taken 
into account. The results are presented in Section[I] We then 
discuss our results in Section and present our conclusions 
in Section [§] 



2 THERMAL COMPTONIZATION IN 
ACCRETING BLACK HOLES 

Table summarizes best-fit plasma parameters of thermal 
Comptonization models fitted to Seyferts and two black-hole 
binaries in the hard state. This model has been found to fit 
the observed spectra rather well. We see that the fitted pa- 
rameters cover a remarkably narrow range of the electron 
temperature, kT c ~ 50-100 keV. A similarly narrow range 
is obtained for the Thomson optical depth; the correspond- 
ing range for spherical geometry is r ~ 2.5-1.5. For clarity, 
the uncertainties are not given in TableQ however, the best- 
fit values from a rather large number of fits all fall within 
the above range, which confirms the statistical significance 
of this result. The points for spherical geometry (and in- 
cluding results of fits of hybrid plasma) are shown in Fig. 
Those points were mostly obtained with the highly accurate 
Comptonization model of Poutanen & Svensson (1996). 

In the case of Seyferts, the data existing so far are 
mostly insufficient to constrain the parameters of the Comp- 
tonizing plasma, with the exception of the Seyfert brightest 
at ~ 100 keV, NGC 4151. To circumvent this problem, av- 
erage spectra from some sets of observations were formed 
(Nandra & Pounds 1994; Zdziarski et al. 1995; Gondek et al. 
1996; Zdziarski, Poutanen, & Johnson 2000, hereafter ZPJ). 
Among them, only ZPJ fit the resulting average spectra with 
thermal Comptonization, and thus we use their results here. 

We note that ZPJ used as free parameters kT e and y 
instead of r. To be able to show the results in our Fig. 
we have refitted the average spectrum of Seyfert Is in the 
kT c -T space with the same assumptions as in ZPJ. Their 
data are from the OSSE detector aboard Compton Gamma 
Ray Observatory, which covers the energy range of above 
50 keV only. Data below 50 keV provide then additional 
constraints. However, the data sets of ZPJ were found to 
be compatible with the range of the X-ray photon spectral 
index of P ~ 1.5-2.3, which corresponds quite well to the 
range observed in Seyferts (e.g., Nandra & Pounds 1994; 



Table 1. Hot plasma parameters in accreting black holes. 



Object 


kT c [keV] 


r a 


Geometry 


u 

Reference 


average Sy 1 


69 


1.6 


sphere 


ZPJ 


average Sy 2 


81 


1.7 


sphere 


ZPJ 


NGC 4151 


73 


1.5 


sphere 


7 




(55 


1.5 


sphere 


ZPJ 




62 


2.0 


hemisphere 


3 


Cyg X-l 


100 


2.0 


hemisphere 


2 




100 


1.3 


sphere 


2 




59 


1.9 


sphere 


1 




58 


2.9 


sphere c 


4 


GX 339-4 


52 


0.9 


slab 


6 




57 


2.0 


sphere 


6 




46 


2.5 


sphere 


5 




43 


2.7 


sphere 


5 




58 


1.9 


sphere 


5 




46 


2.2 


sphere 


5 




76 


1.5 


sphere'" 


5 



a along the radius for sphere/hemispehere, and half-thickness for 
a slab. 

b 1: Frontera et al. (2001); 2: Gierlinski et al. (1997); 3: Johnson et 
al. (1997); 4: McConnell et al. (2002); 5: Wardzinski et al. (2002); 
6: Zdziarski et al. (1998); 7: Zdziarski et al. (2002). 
c A hybrid (Maxwellian with a tail) electron distribution as- 
sumed. 

Zdziarski et al. 1999, 2003). In the contour shown in Fig. Q 
the end with the lowest r and the highest kT c corresponds to 
the softest spectra, with V ~ 2.2-2.3, whereas the opposite 
end corresponds to T ~ 1.5. Thus, the extent of this error 
contour actually corresponds to the range of the observed V. 
It is also remarkable that this error contour agrees relatively 
well with the best-fit parameters of a number of individual 
accreting black holes (for the same geometry). 



3 THE EQUATIONS AND NUMERICAL 
APPROACH 

We concentrate on the inner region of the accretion disc 
since most of the radiation comes from it. We adopt the po- 
tential of Paczynski & Wiita (1980) to mimic the geometry 
of the central Schwarzschild black hole. Steady axisymmet- 
ric and two-temperature assumptions to the accretion flow 
are adopted as usual. A randomly oriented magnetic field 
is assumed to exist in the accretion flow and the magnetic 
pressure is in equilibrium with the gas pressure. We further 
assume that the accretion rate is independent of radius and 
the viscous dissipation mainly heats ions. We will discuss 
the validity of these two assumptions in Section^] 

With these assumptions, we can write down the equa- 
tions describing the dynamics of the hot accretion flow 
(ADAF & LHAF). These equations include the conserva- 
tions of mass and momentum fluxes, and the energy equa- 
tions of ions and electrons. For the case of the standard 
electron energy equation, these equations are the same as 
in Paper I. In particular, cooling due to Comptonization in 
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Figure 1. The observed parameters, electron temperature, kT e 
and Thompson optical depth, t, of a number of accreting black 
holes compared to results of two hot accretion models, ADAF 
and LHAF. The contour corresponds to the average Seyfert- 1 
spectrum of ZPJ. The model points are from Tables [2] and [3] 
but excluding the case of a = 0.01. The solid and dashed lines 
correspond to y = 1 and 0.5, respectively. 



Rout is large, the range is small, hence the effect of OBC is 
not very important. 

We set the T B and T\ at the outer boundary at any values 
for which we can obtain a physical solution. There is no 
difference in this respect between the ADAF and the LHAF. 
Due to the finite range of the allowed OBC, such a solution is 
not unique. However, we confirm that for both ADAFs and 
LHAFs the allowed ranges are narrow, and thus the effect of 
the OBC does not affect our results. The main reason for this 
is that the mass accretion rates we consider are relatively 
high. Another reason, in the case of parameterized energy 
equation, is that the differential equation is replaced by an 
algebraic one. 

We would like to emphasize that the narrowness of 
the range of the OBCs depends on the numerical method 
adopted. Two methods to obtain global solutions of accre- 
tion flows are generally used, namely the shooting and relax- 
ation ones. When the (more accurate) shooting method is 
adopted, as in the present paper, the narrowness of the OBC 
is universal for the global solutions of any accretion models 
(e.g., ADAFs, Nakamura et al. 1997; slim disc, Abramow- 
icz et al. 1988). It does not imply the solution cannot be 
realized in nature. When the OBC of a flow is out of this 
narrow range, it is likely that the flow can adjust itself to 
find the actually existing solution by some other physical 
processes, e.g., thermal conduction, or the solution can be 
weakly time-dependent on the viscous time scale. 



that case is treated using a local energy enhancement fac- 
tor (Dermer, Liang, & Canfield 1991; Esin et al. 1996). On 
the other hand, Comptonization of external soft photons is 
taken into account by using the energy equation of elec- 
trons parametrized by the Compton parameter, analogously 
to Zdziarski (1998). 

The numerical approach is also the same with in Paper 
I. The physical solutions should satisfy the sonic point condi- 
tion, the zero-torque condition at the horizon, and an outer 
boundary condition (OBC). Since we are interested in the 
inner region of the disc, we set the radius of the outer bound- 
ary at R ou t = 10 2 J? S . Our solution is self-consistent since we 
simultaneously solve the coupled radiation-hydrodynamical 
equations including the energy equations of both ions and 
electrons. This ensures that our obtained values of T G and r 
are exact enough to be compared with the observations. 

The OBC is found to play an important role in deter- 
mining some details of the ADAF, such as the value of T e 
(Yuan 1999; Yuan et al. 2000). This is because the differ- 
ential terms such as that describing the advection of energy 
play an important role in the equations, so the description of 
accretion is mathematically a boundary-value problem. Be- 
cause the equations describing ADAFs and LHAFs are com- 
pletely identical, we expect similar effect exists for LHAFs 
as well. We therefore should choose the OBC reasonably 
based on some physical considerations. For example, if the 
hot flow is formed by some processes from an outer cold disc, 
we should consider the transition process. An understanding 
of such process is still lacking. However, it was found that 
physical accretion solutions exist only for a range of T e and 
Ti at the OBC (Yuan 1999). The actual range depends on 
the parameters of the flow. For example, when m is high and 



4 RESULTS 

4.1 Seyfert galaxies 

4- 1.1 The parameterized electron energy equation 

Wc first calculate the critical mass accretion rate, rhi, of an 
ADAF parameterized by the Compton parameter y, assum- 
ing M — 10 8 Mq. The usual advection factor is defined as 
/ = g a dv/9vis- The canonical ADAF has / ~ 1, an ADAF at 
the critical rate has / — 0, while an LHAF has / < 0. Since 
/ is a function of radius, we define the critical rate, rhi, as 
the maximum rate at which / for the corresponding solu- 
tion is still > at all radii, see the dashed curve in Fig.^b). 
With this definition, we find that for R ou t = W 2 R S , rhi of 
ADAFs parameterized by y is 0.04, 0.4, 1, 2, 0.5 for (a,y) = 
(0.01,1), (0.05,1), (0.1,1), (0.3,1), (0.1,0.5), respectively. 
Very approximately, rhi — y(a/0.1) 1A for a<0.1, similar 
to rhi ~ lOy - 6 ^ 1 ' 4 obtained by Zdziarski (1998). 

In Paper I, we found that LHAF has two possible struc- 
tures depending on the accretion rate. When rhi < rh < rhi, 
the accretion flow remains hot throughout the disc. Here 
m2 ~ 5mi , and its exact value depends on the flow parame- 
ters and the OBC. When rfi2 < m < 10, the Coulomb cooling 
of the ions becomes so efficient (due to the high density) 
within certain radius that the accretion flow collapses onto 
the equatorial plane and forms a cold, optically-thick, annu- 
lus. Similarly in the present case of the parameterized energy 
equation, there are also two types of LHAF, one remaining 
hot throughout the disc and the other collapsing and form- 
ing a cold optically-thick annulus within a certain radius. 
Note that we assume that the cooling of the hot plasma by 
the optically-thick emission of the cold annulus is included 
in our parametrized electron energy equation. 
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Figure 2. The radial variation of the radiation rate, q le Iir (in dimensionless units with M = G = c = 1 and multiplied by a factor of 
10 17 ), the advection factor, <j a dv/<7vis> T e and t, for two hot disc solutions with the parameterized electron energy equation. The dashed 
curves correspond to an ADAF solution, with a = 0.3, y = 0.89, m = 2 (f» mi), and the OBC at T- = 8 X 10 9 K, T c = 4 X 10 9 K. The 
solid curves correspond to an LHAF solution, with a = 0.3, y = 0.89, rh = 4, and the OBC at T; = 9 X 10 9 K, T c = 2.8 X 10 9 K. 



Fig- |U shows two examples of the global solution, with 
the dashed curve denoting a critical ADAF and the solid 
curve an LHAF. The four panels show the radial profiles of 
T c , q\ c r 2 H, r, and q a dv/q v i B - We see that the T e and r de- 
pend on radii. This presents us an issue of which values of 
T c and r are most representative for the emitted spectra, 
and can then be compared with results of fits to data with 
Comptonization in a uniform plasma. We see in Fig. [5Ja) 
that the emission rate per logarithm of radius, proportional 
to qi e r 2 H, show rather strong maxima, both for the ADAF 
and the LHAF. We find it to be generally the case, both for 
the standard and the parameterized electron energy equa- 
tion. Therefore, we compare the values of T c and r at this 
maximum with the observational data. 

Since the maximum peak is generally sharp, a uniform 
slab appears to be an approximation to the flow geometry 
somewhat worse than that of a sphere. On the other hand, 
the actual flow geometry is still flattened, and the values of 
t obtained in the spherical geometry should be somewhat 
decreased when compared to the hot-flow results. In the case 
of a uniform slab, the fitted r of the half-thickness is lower 
by a factor of > 2 from the radial r fitted in the spherical 
geometry (e.g., ZPJ). We thus estimate that the values of t 
obtained from fits in spherical geometry should be reduced 
by a factor of ~ 1.5 when compared with the r of a hot flow 
calculated as above. Note that this correction also implies 
that the values of y used in the hot flow models should be 
correspondingly reduced to correspond to the data. 

Our results, including the model luminosity, for the pa- 
rameterized electron energy equation are given in Table |5J 
and the obtained values of kT B and r are also shown in Fig. 
The luminosity is calculated assuming the electrons radi- 
ate away all the energy obtained from Coulomb collisions, 



which is a good approximation at high accretion rates. The 
luminosity for solutions with a cold annulus is given only 
for the hot flow, i.e., not including that of the annulus. We 
consider a = 0.01, 0.1, 0.3, and y = 0.5, 0.89, 1. The middle 
value of y corresponds to the best fit to the average Seyfert-1 
spectrum by ZPJ (but without correcting for the difference 
between geometries of the flow geometry and of a sphere). 
The considered range of y roughly corresponds to the range 
of the 2-10 keV spectral index most common in Seyferts, 
T ~ 1.7-2.1 (Nandra & Pounds 1994; Zdziarski et al. 1999). 

Our goal is to check which accretion models yield T c 
and t in agreement with the observaional data for luminous 
sources. For the ADAF, we only show solutions with m ~ 
rhi since solutions with a lower m would yield values of T e 
clearly much higher than those observed. 

For a = 0.1, 0.3 and y = 0.5, 0.89, 1, we find that T c 
for parameterized ADAFs is > 10 9 ' 3 K, while it is 10 8 ' 9 K < 
T c < 10 9 ' 2 K for parameterized LHAFs. The agreement be- 
tween the predictions of the LHAF solution and the ob- 
servations is very good. We therefore conclude that these 
LHAF solutions (with additional cooling responsible for the 
observed values of y) are very likely to correspond to the 
actual accretion flows in Seyfert-1 galaxies. 

We have also considered a — 0.01. We find that in this 
case neither ADAF nor LHAF can be reconciled with the 
data because of the too high T e predicted. This provides an 
interesting observational constraint on the viscous parame- 
ter, a > 0.1. If the viscous stress is caused by magnetic field 
via the magneto-rotational instability, the numerical simu- 
lation by Hawley, Gammie & Balbus (1996) imply a S> 0.01 
(with a ~ 0.4 for equipartition magnetic field), in agreement 
with our constraint. 

In our calculations, we assume that the Compton pa- 
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Table 2. Model parameters of some parameterized hot disc solutions. 



Number 


Type a 


a 


V 


TO 


log(T c [K]) 


T 


L/L E [%] 


1 


A Ft A F 1 


U.o 


u.oy 


1(< mi) 


O 

y .04 


U.oo 


U.o 


O 

z 


A Ft A F 


U.o 


U.oy 




y . o 


U.OO 


o . 




r IT AT? 


U.o 


n so 
U.oy 


A 
1 


n i t; 

y . io 


n oq 
U.yo 


D. 1 


A 
1 


T TT A T? 


U.o 


U.oy 


1 


y.ui 


1 OQ 


i n o 


r 




T ITAT? 


U.o 


U.oy 


C 
O 


y .U4 


1 1 o 

1. ly 


R p. 
U . u 


6 


LHAF 


0.3 


0.89 


5 


9.03 


1.23 


5.0 


V 


LHAF 


0.3 


0.89 


6 


9.07 


1.13 


4.0 


8 


LHAF 


0.3 


0.89 


6 


9.06 


1.14 


3.5 


9 


LHAF b 


0.3 


0.89 


5 


9.074 


1.11 


2.5 


10 


LHAF 6 


0.3 


0.89 


6 


9.069 


1.13 


3.2 


11 


ADAF 


0.3 


0.5 


1(< mi) 


9.3 


0.37 


1.0 


12 


ADAF 


0.3 


0.5 


1(< mi) 


9.4 


0.29 


1.0 


13 


LHAF 


0.3 


0.5 


2(> mi) 


8.95 


0.84 


8.3 


11 


LHAF 


0.3 


0.5 


4 


8.96 


0.82 


3.0 


15 


LHAF 


0.3 


0.5 


5 


8.99 


0.76 


2.5 


16 


LHAF 


0.3 


0.5 


8 


9.0 


0.74 


3.2 


17 


ADAF 


0.3 


1.0 


1(< mi) 


9.6 


0.37 


0.7 


18 


ADAF 


0.3 


1.0 


2 (ss mi) 


9.36 


0.65 


3.5 


19 


LHAF 


0.3 


1.0 


4 


9.03 


1.38 


6.74 


20 


LHAF 


0.3 


1.0 


4 


9.13 


1.1 


11.1 


21 


LHAF 


0.3 


1.0 


4 


9.02 


1.4 


7.35 


22 


LHAF 


0.3 


1.0 


5 


9.05 


1.34 


7.8 


23 


ADAF 


0.1 


0.89 


0.8(s3 mi) 


9.4 


0.52 


0.9 


24 


LHAF 


0.1 


0.89 


2 


9.1 


1.03 


6.3 


25 


LHAF 


0.1 


0.89 


2 


9.4 


0.52 


1.6 


26 


LHAF 


0.1 


0.89 


4 


9.16 


0.92 


8.2 


27 


LHAF 6 


0.1 


0.89 


3 


9.03 


1.23 


1.6 


28 


LHAF 6 


0.1 


0.89 


4 


9.07 


1.13 


2.7 


29 


LHAF 


0.01 


0.5 


0.1 


9.6 


0.19 


0.05 


30 


LHAF 


0.01 


1.0 


0.1 


9.6 


0.37 


0.19 



a Models with the same a, y, and rh have different outer boundary condition. 
^ There is a cold annulus transition within a certain radius in these LHAF models. 
In these cases, L is given for the emission of the hot flow only. 



rameter is independent of the radius. While this is a simplifi- 
cation, we do not expect it to substantially affect our results. 
This is because in the case of a variable y, we can simply 
set the observational value of y (e.g., « 0.89) at the radius 
where most of radiation comes from. Then, since m oc rv r r, 
the only way that a variable y could affect our results is 
through modifying the radial velocity, v T , which is not signif- 
icant. The main physical reason why LHAF solutions yield 
lower T e compared to the ADAF is because the former cor- 
responds to higher accretion rates and thus higher r. 



4-1.2 The standard electron energy equation 

As discussed in Section 1, global solutions including emis- 
sion of soft seed photons by all possible sources is both diffi- 
cult and underdetermined given our present understanding 
of the flow dynamics. Therefore, we have parametrized the 
energy equation in Section 14.1.11 by the observational con- 
straint given by the X-ray slope. On the other hand, it is also 
important to check how much the presence of the additional 
soft photons affects our solutions. Thus, here we calculate T c 



and t of the ADAF/LHAF solutions with a standard energy 
equation in which no external soft photons are included in 
the Comptonization process. 

The results are shown in Table |S| Again we find that 
the obtained values of T c in the ADAF case are too high, 
and thus can be ruled out for luminous sources. We also find 
that the LHAF without a cold annulus provides marginally 
appropriate values of T e and r in some cases, and not appro- 
priate in some other cases (e.g., model 9 in Table El . This 
indicates the importance of the Comptonization of external 
soft photons. On the other hand, LHAFs with a transition 
to a cold annulus always give appropriate values of T c and 
t. In addition, from the calculated T c and r of both ADAF 
and LHAF, we calculated y and find that it is generally < 1. 
Note that we do not include here the effect of cooling by the 
cold annulus. 

We also present our results in Fig. Comparing with 
the contour which corresponds to the average Seyfert-1 spec- 
trum (ZPJ), we find that the LHAF is much more favoured 
than the ADAF as the accretion disc model describing the 
X-ray emission of Seyferts. 
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Table 3. Model parameters of some standard hot disc solutions. 



Model 


Model type a 


a 


lit 


l°g(T c [K]) 


T 


y b 


L/Le [%} 


1 


AiJAr 


0.3 


1(R2 mi) 


9.44 


0.38 


0.71 


0.9 


2 


ADAF 


0.3 


1 


9.46 


0.34 


0.66 


0.7 


3 


ADAF 


0.3 


1 


9.29 


0.7 


0.92 


3.2 


1 


ADAF 


0.3 


1 


9.26 


0.73 


0.90 


3.9 


5 


LHAF 


0.3 


2 


9.22 


0.81 


0.91 


5.9 


6 


LHAF 


0.3 


.3 


9.24 


0.76 


0.89 


6.0 


7 


ADAF 


0.1 


0.1(< mi) 


9.85 


0.05 


0.24 


0.004 


8 


ADAF 


0.1 


0.5(« mi) 


9.5 


0.3 


0.64 


0.36 


9 


LHAF 


0.1 


1 


9.5 


0.2 


0.43 


0.2 


10 


LHAF C 


0.1 


3 


8.9 


1.35 


0.72 


4.4 


11 


LHAF C 


0.1 


5 


9.08 


1.12 


0.91 


2.4 


12 


LHAF C 


0.1 


10 


8.9 


0.83 


0.45 


1.3 



a Models with the same a, and m have different outer boundary condition. 
^ Calculated from the T c and r (i.e., not a free parameter). 

c There is a cold annulus transition within a certain radius in these LHAF models. In these 
cases, L is given for the emission of the hot flow only. 



4.2 The hard state of X-ray binaries 

In Section 14. II we have applied our model to AGNs. Here, 
we extend our calculations to the case of stellar-mass black 
holes, assuming M — IOMq. We find that the results of 
our models are virtually independent of the mass. On the 
other hand, the data for black-hole binaries are usually much 
better than those for Seyferts. For example, the values of T e 
are determined rather accurately for some hard states of 
Cyg X-l and GX 339-4, with T c w 50-60 keV (Frontera et 
al. 2001; Zdziarski et al. 1998; Wardzihski et al. 2002). Such 
low values of T e can hardly be reached in the ADAF case, 
while they are within the range of T e predicted by the LHAF 
model (Table (A specific application of the LHAF model 
to an X-ray binary, XTE J1118+480, has recently been given 
by Yuan, Cui & Narayan 2004.) 

In addition to the value of T e , another very strong ar- 
gument in favour of the LHAF is the very large bolometric 
luminosity of transient (with a low-mass companion) black- 
hole binaries in the hard state. During the rising phase of 
an outburst, luminosities as high as ~ 0.2Le are observed 
(e.g., in XTE J1550-564, Done & Gierlihski 2003; GX 339- 
4, Zdziarski et al. 2004, see also Nowak 1995; Maccarone 
2003). Such luminosities cannot be achieved in the ADAF 
model, but they can (within a factor of two) be obtained by 
the LHAF model (see Table 0. 

Narayan (1996) proposed that the variety of spectral 
states of black-hole X-ray binaries can be understood as a 
sequence of thin disc plus ADAF models with varying rh 
and the transition radius. This idea was developed by Esin, 
McClintock & Narayan (1997) and Esin et al. (1998). In 
their work, the hard state is associated with ADAFs with 
m<rhi. The advection factor /, which is a function of ra- 
dius, is instead set to a value averaged over the whole ac- 
cretion flow. Then / is found to have a low but still posi- 
tive value, / ~ 0.3, at m ~ rhi. The model typically pre- 
dicts T c > 10 9 K for accretion flows within ~ 100i? s , and 
T e > 10 9 ' 1 K for r < 10i? s where most of the radiation comes 
from (see Fig. 3b in Esin et al. 1997). This is in good agree- 



ment with our results, in which we also find that an ADAF 
predicts T e > 10 9 ' 1 K. While such an ADAF model can ex- 
plain the spectra of the hard state of some X-ray binaries 
very well, it is challenging to fit other spectra with lower en- 
ergy cutoff such as Cyg X-l and GX 339-4 because a lower 
T e < 10 9 K is required there. 



5 DISCUSSION 

A caveat for our results concerns the use of the one-zone ap- 
proximation of Comptonization in comparisons with predic- 
tion for T c and r of accretion flow models. We use it because 
most of the available Comptonization fits to data in astro- 
physical literature use the one-zone approach (i.e., a uniform 
electron temperature in a simple geometry). Thus, in order 
to compare our results with the data, we neglect the radial 
dependences of the temperature and density and assume 
that most of the radiation comes from a single radius. We 
intend to extend our calculation in the future by including 
the profiles of the electron temperature and density in the 
accretion flow and the non-local nature of Comptonization. 
On the other hand, the above approximation is not used in 
comparisons of the predicted and observed luminosities of 
black-hole binaries. 

In our calculations, we assume that the accretion rate 
is independent of radius. Over the past few years, hydrody- 
namic and magnetohydrodynamic simulations (e.g., Stone, 
Pringle & Begelman 1999; Hawley & Balbus 2002; Igumen- 
shchev et al. 2003) and analytical work (Narayan & Yi 1994; 
Blandford & Begelman 1999; Narayan et al. 2000; Quataert 

6 Gruzinov 2000) indicate that only a fraction of the gas 
that is available at large radius in the accretion flow may 
actually accrete onto the black hole. The rest of the gas is 
either ejected from the flow or is prevented from accreting 
by convective motions. The former is due to the positive 
sign of the Bernoulli parameter of the accretion flow while 
the latter is due to the the accretion flow being convection- 
unstable. However, all the above results are for very low 
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accretion rates. When the accretion rate is high, rh > rh\ , as 
in our case, the Bernoulli parameter is negative in general, 
and the flow is likely to be convection-stable since the en- 
tropy of an LHAF decreases rather than increases towards 
the smaller radii (Paper I) . Therefore, the constant accretion 
rate is likely a correct assumption in our case. 

We also assume that the viscous dissipation mainly heat 
ions. Depending on some unknown details of microphysics, 
it is possible that a large fraction of the viscous dissipation 
also heat electrons (Quataert & Gruzinov 1999). If this is 
the case, one effect would be that the value of T e predicted 
in both ADAF and LHAF will be somewhat larger for a 
fixed rh. This will make the ADAF even worse for describing 
the X-ray emission of Seyfert-1 galaxies and the hard state 
of X-ray binaries. Another effect of electron heating would 
be that the value of the critical accretion rate of ADAF, 
mi, becomes smaller due to the weaker viscous dissipation 
heating of ions. In this case, the lowest T c an ADAF can 
produce will become higher, which again makes the ADAF 
worse as the model of Seyfert-1 galaxies and X-ray binaries 
in the hard state and implies the occurrence of LHAFs. 

There are two possible origins for the hot accretion 
flows. One is through the transition from an outer standard 
thin disc. Several mechanisms have been suggested for the 
transition from the outer cold disc to the inner hot disc. One 
is the evaporation of the cold disc due to thermal conduc- 
tion (Meyer & Meyer-Hofmeister 1994; Meyer, Liu, & Meyer- 
Hofmeister 2000; Rozariska & Czerny 2000). The second is 
turbulent diffusive heat transport (Honma 1996; Manmoto 
et al. 2000; Manmoto & Kato 2000). The third one is the sec- 
ular instability present in the radiation pressure-dominated 
inner region of the standard thin disc (Lightman & Eardley 
1974; Gammie 1998) . If the accretion flow is one-phase, then 
the accretion mode within R ou t is simply determined by the 
accretion rate, rh. If rh < rhi, it is an ADAF, otherwise it 
is an LHAF. If the transition is due to secular instability, 
the instability may result in the formation of cold clumps 
embedded in the hot flow. In this case, the accretion rate 
and accretion mode in the hot phase may depend on the 
interchange of matter and energy between the cold and hot 
phases. The emission from the cold clumps will supply ad- 
ditional soft photons as the seeds for Comptonization. 

The second possible origin for hot accretion flow is that 
the accretion flow is already hot at large radii (e.g., Shlos- 
man, Begelman & Frank 1990). First, there exists such a hot 
branch of accretion solution even for the accretion rates as 
high as Eddington (Paper I). Second, in the case of AGNs, 
there is strong observational evidence that hot cooling flows 
carry large amount of hot gas into the central region of 
AGNs (Sarazin 1986; Arnaud 1988) which could serve as the 
accretion material. The hot ISM is also a source of accret- 
ing material. In the case of X-ray binaries, if the accretion 
comes from the supersonic stellar wind from the companion 
star, the accretion flow may start out hot due to the heating 
of a bow shock. 



6 CONCLUSIONS 

In this paper, we have investigated the hot accretion disc 
model applied to the X-ray emission of Seyfert galaxies and 
black-hole binaries in the hard state. Especially, we con- 



sider the luminous hot accretion flow recently found (Paper 
I) — LHAF — and compare it with the ADAF. We numeri- 
cally solve the radiation hydrodynamical accretion equation 
to obtain the electron temperature and the Thomson op- 
tical depth of these two models at the radius where most 
of the emission comes from. We adopt two forms of the 
electron energy equation, with one being the standard in 
the sense that only synchrotron, bremsstrahlung and their 
Comptonization are included as the emission mechanisms, 
with the other parameterized by the Compton parameter, y, 
to approximate the actual cooling process of electrons, e.g., 
external soft seed photon and feedback process between cold 
and hot components. 

Comparing our results to those obtained by fitting the 
observed average spectra of Seyfert-1 galaxies and black- 
hole binaries by thermal Comptonization models, we find: 
(i) the ADAF with the standard electron energy equation 
is ruled out due to its T c being too high; (ii) the ADAF 
with external soft photon as seed photon of Comptonization 
is also not likely to be responsible for the X-rays; (iii) the 
LHAF with standard electron energy equation is marginally 
feasible in the sense that T c is in the edge of the observed 
parameter space; (iv) the most possible accretion disc model 
for the X-ray emission of Seyfert-1 galaxies is an LHAF with 
the electron energy equation including additional sources of 
seed photons (parametrized by y); (v) the high bolometric 
luminosity of some black hole X-ray binaries in the hard 
states can be achieved in an LHAF but not an ADAF. 
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